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ABSTRACT 

Luminous X-ray gas coronae in the dark matter halos of massive spiral galaxies are a fundamental 
prediction of structure formation models, yet such coronae remained essentially unexplored. In this 
paper, for the very first time, we detect and characterize extended hot X-ray coronae beyond the 
optical disks of two "normal" massive spiral galaxies, NGC1961 and NGC6753. Based on XMM- 
Newton X-ray observations, we detect hot gaseous emission extending out to ~60 kpc around both 
galaxies - well beyond their optical radii. The hot gas, whose best-fit temperature is kT ~ 0.6 keV and 
abundance is ~0.1 Solar, appears to have a fairly uniform distribution, hinting that the quasi-static 
gas resides in hydrostatic equilibrium in the potential well of the galaxies. The bolometric luminosity 
of the hot gas in the (0.05 — 0.15)r2oo region, where r2oo is the virial radius, is x 10*° erg s~^ 
for both NGC1961 and NGC6753. We derive the baryon mass fractions of NGC1961 and NGC6753, 
and obtain /b,NGCi96i ~ 0.11 and /b,NGC6753 ^ 0.09, which values fall short of the cosmic baryon 
fraction. The detected hot X-ray coronae around NGC1961 and NGC6753 offer an excellent basis to 
probe structure formation simulations. To this end, the observations are confronted with the recently 
developed moving mesh code AREPO and the traditionally used smoothed particle hydrodynamics 
code GADGET. The implemented subresolution physics and the gravity solver arc fully identical in the 
two codes, but they use different methods to solve the hydrodynamical equations. We conclude that, 
while neither model gives a perfect description, the observed luminosities, gas masses, and abundances 
favor the AREPO code. Moreover the shape of the observed density profiles are also well reproduced 
by the moving mesh code within ~0.4r2oo- However, neither model incorporates efficient feedback 
from supermassive black holes or supernovae, which could alter the simulated properties of the X-ray 
coronae. With the further advance of numerical models, the present observations will play a crucial 
role in constraining the feedback effects in structure formation simulations. 

Subject headings: galaxies: individual (NGC1961, NGC6753) — galaxies: spiral — galaxies: ISM — 
X-rays: galaxies — X-rays: general — X-rays: ISM 



1. INTRODUCTION 

The presence of hot gaseous coronae in dark matter 
halos of massive galaxies is a fundamenta l prediction of 
structure formation models (jWhite fc Ree s 1978). This 
gas collapse scenario in a cold dark matter dominated 
Universe was applied bv I White fc Reed ()199H ). who pre- 
dicted that dark matter halos - surrounding all types of 
present epoch massive galaxies - should have an associ- 
ated X-ray luminous corona. 

X-ray luminous, hot gaseous coronae are nearly 
ubiquitous for massive early- type galaxies. First 
discovered in Ein stein Observatory observations 
(jForman et al.l 119851 ). these corona e have been exten- 
sively studied fe.g. ISarazin et all 120011: IJones et al 
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200a [Mathews fc Brighentil 120031: IBogdan fc Gilfanov 
2011b[ ). Early-type galaxies exhibit at least three X-ray 
emitting components. (1) A well-studied population 
of low-mass X-ray binaries (LMXBs), which add a 
notable contribution to the total X-ray emission of 
galaxies (|Gilfanovll200"l IZhang et al.llMl) : (2) emission 
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from active binaries (ABs) and cataclysmic variables 
(CVs), whose total emission is proportional to the 
stellar mass of the galaxy (e.g. iRevnivtsev et al.l 120081 : 
IBogdan fc Gilfanovl |20TTH ). and (3) thermal, sub-keV 
gas, produced from the galaxys own constituent stars or 
from infall of group/cluster gas in central bright cluster 
galaxies. For massive early-type galaxies, it is difficult 
to detect the more extended atmospheres of hot gas, 
since they either have bright coronae from stellar mass 
loss and/or lie in larger dark matter halos - groups and 
clusters - that have their own gaseous atmospheres. 
Thus, elliptical galaxies are not suitable for probing 
the luminous X-ray gas predicted around individual 
galaxies. 

However, late-type, spiral galaxies, also are predicted 
to have luminous X-ray coronae and c an test galaxy for- 
mation scenarios (|White fc Reeslll991| ). In fact, the most 
promising probe of structure formation models is via 
the study of spiral galaxies. The observed X-ray emis- 
sion from spiral galaxies also consists of multipl e com- 
ponents: (1) In addition to LMXBs ()Gilfanovl [20041. 
bright high- mass X-ray binaries (HMXBs), associated 
with star-formation, are located in the star-forming re- 
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TABLE 1 

The properties of the sample galaxies. 

Name Distance 1' scale Njj Lk Mi,/LK M* SFR Morph. r2oo Tobs 

(Mpc) (kpc) (cm-2) (Lk.©) (Mq/Lk.q) (M©) (Mq yr"!) type (kpc) (ks) 

(1) (2) (3) (4) {h) {&) {]) (8) (9) (10) 

NGC1961 55.8 16.23 8.4 x lO^o 5.4 x lO" 0.78 4.2 x 10" 15.5 SAB(rs)c 470 73.7 

NGC6753 43.6 12.67 5.5 x 10^0 3.9 x lO" 0.81 3.2 x lO" 11.8 (R)SA(r)b 440 73.9 



Note. Columns arc as follows. (1) Distance taken from NED |htt p : / / ned ww w .ipac.caltech.edu/ ). (2) 1' scale at the applied dista nce. (3) 
Galactic absorption (Dickey & Lockman 1990). (4) Total K-ban d luminosity. (5j K-band mass-to-light ratios computed from B ell et al.l 
l |2003 ) using the B ~ V color indices of galaxies llde Vaucou leurs et al.. ,1991il . (6) Total stellar mass based on the K-band luminosity and 
the K-band mass-to-light ratios. (7) Star formation rate computed from the IRAS 60 /xm and 100 ^m flux, details are given in Section 
13.31 (8) Morphological type, taken from NED. (9) Virial radius of the galaxies, estimated from the maximum rotation velocity. (10) Total 
XMM-Newton exposure time. 



gions of spiral galaxies ()Grimin "eFall 12001 : (2) Also 
an important contributor to the X-ray flux of spiral 
galaxies is the collective emission of faint compact X- 
ray sources, which contribute to the unresolved, diffuse 
component. The populat ion of ABs and CVs fol low the 
stellar light distribution (jRevnivtsev et al.ll2008[ ). while 
the emission from young stars and young stellar objects 
is proport ional with the star-format ion rates of the host 
galaxy (Bogd an fc Gilfanovl [20TIbl) . (3) Additionally, 
spiral galaxies host at least moder ate amounts of ion- 
ized ga s with sub-keV temperatures (jBogdan fc Gilfanovl 
l2011bf ). These components, enumerated above, provide 
the bulk of the observed emission within the optical ex- 
tent of spiral galaxies and are very difhcult to model ac- 
curately as an integrated component to detect the faint 
outer corona projected onto this complex region. There- 
fore, to probe the structure formation models, the hot 
X-ray coronae must be explored beyond the optical ex- 
tent of spiral galaxies. 

A major advantage of spirals is their location: whereas 
massive ellipticals lie either in the center of galaxy groups 
and clusters and/or in rich environments, it is possible to 
find relatively isolated massive spiral galaxies. Moreover, 
to study the X-ray coronae around spiral galaxies, a qui- 
escent environment is essential. Hence, starburst galax- 
ies and galaxies undergoing mergers are not suitable to 
characterize the extraplanar X-ray emission. Thanks to 
UV, far-infrared, and Ha imaging, large star-formation 
rates, hence starburst galaxies, can be identified, and 
merging galaxies can be recognized based on their op- 
tical appearance. Thus, isolated massive spiral galaxies 
with moderately low star formation rates and relatively 
undisturbed morphology are the ideal targets for explor- 
ing the extraplanar X-ray coronae, thereby testing struc- 
ture formation models. 

Our current theoretical view about hot gaseous coro- 
nae around massive spirals is based on numerical simula- 
tions of galaxy formation. The first analytic analysis of 
galaxy forma tion in the cold dark m atter cosmogony was 
presented bv iWhite fc ReesI ()1991[ ). who suggested that 
X-ray coronae around massive galaxies should be ubiqui- 
tous. For example, their calculations argued that a halo, 
characterized by 300 km s''^ circular velocity, would have 
an X-ray luminosity of 3 x 10"*^ erg s""'^ and a gas tem- 
perature of about 0.3 keV . How ever, the luminosities 
predicted bv IWhite fc ReesI (jl991l ) are too high, because 
they neglected gas ejection and assumed that the hot gas 
follows the dark matter distribution . The hydrodynam- 
ical simulations of iToft et al.l (|2002D predicted less mas- 



sive and less luminous X-ray coronae than the earlier an- 
alytic models. In particular, the simulated X-ray coronae 
are about two orders of magnitude fainter than the orig- 
inal IWhite fc ReesI ([l99l[ ) predictions. Th e lower pre- 
dicted luminosities in the iToft et all (|2002[ ) simulations 
are due to the absence of efficient feedback, which re- 
sulted in very massive stellar components, consequently 
much les s massive and less l uminous X-ray coronae. Most 
recently, iCrain et al.l (|2010[ ) studied the properties of ro- 
tationally supported systems with an extensive simula- 
tion of disk-like galaxies. The simulation, that was car- 
ried out in a large cosmological volume, incorporates ef- 
ficient feedback from supernovae that prevented the con- 
version of halo gas into stars. The supernova feedback 
also altered the hot gas entropy profile by ejecting a frac- 
tion of the baryons to larger radii, and by preferentially 
ejecting the lowest entropy gas to suppress the central 
density profile. 

The primary goals of the present paper are twofold. 
First, we aim to detect and characterize the hot X-ray 
coronae beyond the optical disks of two massive spiral 
galaxies, NGC1961 and NGC6753, based on moderately 
deep XMM-Newton X-ray observations. Second, we aim 
to confront the observed physical properties of the X- 
ray coronae with those predicted by modern structure 
formation models. 

The paper is structured as follows. In Section 2 we 
describe the previously available sample of spiral galax- 
ies and introduce the presently studied NGC1961 and 
NGC6753. In Section 3 we describe the reduction of the 
X-ray and infrared data. The results, namely images, 
surface brightness profiles, and spectra, are presented in 
Section 4. The observed characteristics of the X-ray coro- 
nae are confronted with structure formation simulations 
in Section 5. We discuss our results in Section 6. We 
summarize in Section 7. 

2. PREVIOUS OBSERVATIONS AND 
THE PRESENT SAMPLE OF SPIRAL GALAXIES 

2.1. Previous observations of spiral galaxies 

Prior to the XMM-Newton observations of NGC1961 
and NGC6753, the data available for studying the ex- 
tended X-ray coronae of spiral galaxies was very scarce 
for the following four reasons. First, most studied spi- 
ral galaxies have stellar masses of ^ 3 x 10^^ Mq, 
implying rather faint X-ray coronae beyond their opti- 
cal extent (Section [5.21) . Therefore, the predicted X-ray 
coronae are not detected: only upper limits are obtained, 
which means that the gaseous emission cannot be char- 
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Fig. 1. — Left: 0.3 — 2 keV band raw XMM-Newton image of a 5' X 5' region (81 X 81 kpc) around NGC1961. The net exposure time 
of the galaxy after flare filtering is 31.0, 34.1, and 21.3 ks for MOSl, MOS2, and PN, respectively. The bright X-ray sources beyond the 
optical radius are cosmic X-ray background sources. R ight: DSS R-band image of the same region. The somewhat disturbed morphology 
may be the result of a recent minor merger ICombes et al.ll2009f) . 



acterized. Examples for such spiral galaxies are those 
studied bv lBenson et"all (|2000t ) based on ROSAT PSPC 
data, or the galaxies investigated by iRasmussen et al.l 
(|2009f) based on Chandra observations. Second, most 
galaxies with detected X-ray coronae are not normal spi- 
ral galaxies, but they are either starburst/post starburst 
systems, active galactic nuclei (AGN), or spirals whose 
cold interstellar media have been heated and sometimes 
stripped by the interaction with the intracluster gas. 
Therefore, the starburst galax ies, including the famou s 
M82 and NGC253, studied bv I Strickland et al.l (l2004ah . 
or the galaxies lo cated in rich cluster environment in the 
ISun et all ()2007D sample cannot be used to probe struc- 
ture formation models. Third, extra-planar emission is 
detected in a number of normal, edge-on spiral galaxies, 
which X-ray emission, however, is not extended to large 
radii but lies close to the disk (e.g.|Wang 2005). The ob- 
served extraplanar X-ray emission in these spirals is not 
likely that arising from an extended corona of infalling 
gas, but it is mos t likely a sup e rnova (SN) driven outflow, 
as discussed by iDavid et alj ()2006[ ) for low luminosity 
early- type galaxies and as observed in the bulge of M31 
(|Li fc Wandl2007tlBogdan fc Gilfanovll2008l) . Indeed, de- 
tailed studies of the extraplanar emission demonstrated 
that it is confi ned to within 10 kpc ar ound NGC577 5 
(|Li et al.ll2008l ) and 23 kpc around M104 (|Li et al.ll20lH ). 
and can be interp r eted a s gas driven up from the disk. 
Finally, iDai et ahl (|2011[ ) reported the presence of hot 
gaseous corona around UGC12591, a rapidly rotating 
SO/a (NED) galaxy at 100 Mpc distance. However, 
UGC 12591 is not a definite test of structure formation 
models since 1) it hosts a very luminous AGN, 2) the 
X-ray data is not sufficient to explore the outer regions 
of the corona, hence only the X-ray emission associated 



with the stellar body can be studied, and 3) IDai et ahl 
(|2011f ) employs an overly simplistic spectral fit, which 
does not properly account for various X-ray emitting 
components associated with the galaxy. 

A massive spiral galaxy of special note is NGC1961, 
which is also discussed in the present paper. NGC1961 
was recently observ ed by Chandra for 130 ks 
(jAnderson et al.|[201lD . Despite the deep total exposure, 
the four separate pointings had different aim points and 
roll angles, aiming to study the X-ray gas at very large 
radii from the galaxy. Hence, the exposure for the hot X- 
ray corona, just beyond the optical radii, is very limited. 
Although the data indicates the presence of an extended 
X-ray corona out to 40 — 50 kpc, the Chandra data does 
not allow the characterizatio n of the X-ray corona . In 
fact, the spectral analysis of lAnderson et ahl ()201lD in- 
cluded the bright X-ray emission from the galaxy itself, 
which was then then modeled to extract the spectrum 
of the diffuse corona. Consequently, their derived results 
rely on an assumed abundance and temperature, and the 
extrapolated gas mass is highly uncertain. 

2.2. Spiral galaxies analyzed in this paper 

To detect and characterize hot X-ray coronae with 
present-day X-ray telescopes, the targeted galaxies must 
fulfill the following three criteria. First, the selected 
spiral galaxies must be massive. More massive galax- 
ies are predicted to have more luminous X-ray coronae, 
thereby facilitating detections with high signal-to-noise 
ratio. Second, relatively relaxed galaxies are required in 
a clean environment, that is they should not be star- 
burst systems, should not be undergoing an interaction 
or merger, should not include a signigicant AGN, and 
should not be located in a rich galaxy cluster. Finally, 
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Fig. 2. — Left: 0.3 - 2 kcV band raw XMM-Newton image of a 5' X 5' region (63 X 63 kpc) around NGC6753. The net exposure time 
of the galaxy after flare filtering is 58.7, 59.2, and 27.9 ks for MOSl, MOS2, and PN, respectively. The bright X-ray sources beyond the 
optical radius are cosmic X-ray background sources. Right: DSS R-band image of the same region. The symmetrical morphology of the 
galaxy indicates its undisturbed nature. 



TABLE 2 

The list of analyzed XMM-Newton observations. 



Galaxy 


XMM-Newton 


Centre 


MOSl clean 


MOS2 clean 


PN clean 


Observation 




Obs ID 


coordinates 


exposure time (ks) 


exposure time (ks) 


exposure time (ks) 


Date 


NGC1961 


0673170101 


05h42m04.60s -|-69d22'42.0" 


16.0 


19.1 


12.7 


2011 Aug 31 


NGC1961 


0673170301 


05h42m04.60s -f 69d22'42.0" 


15.0 


15.0 


8.6 


2011 Sep 14 


NGC6753 


0673170201 


19hllm23.60s -57d02'58.0" 


58.7 


59.2 


27.9 


2012 Apr 21 



the radial range of about 25 — 100 kpc should fit in the 
XMM-Newton field-of-view (FOV). The selected galax- 
ies, NGC1961 and NGC6753 fulfill all above listed cri- 
teria. They are fairly massive (M^, > 3 x 10^^ M0), 
are quiescent, and are located in poor galaxy groups 
(jFouque et al.lll99"^ . and are sufficiently distant. Thus, 
NGC196f and NGC6753 are excellent targets to test 
galaxy formation scenarios. The major physical prop- 
erties of the studied galaxies are listed in Table [TJ 

To estimate the virial mass (M200) of the sample galax- 
ies, we use the baryonic Tully-Fisher relation for th e cold 
dark matter cosmogony, M200 k Vj^fx (Navarr o et al.l 
119971) . where Vmax is the maximum rotational velocity 
of the galaxies. The Knax values are taken from the 
HyperLedcQ catalog (jPaturel et al.l [2003) , which gives 
V^^y, = 447.9 ± 14.6 km s'^ for NGC1961 and Fmax = 
395.0 ± 21.1 km s'^ for NGC6753 Usiiig results of nu- 
merical simulations ()Navarro et al.l [l997 ). we estimate 
the virial mass of ~1.2 x 10^^ Mq and ~1.0 x 10^^ Mq for 
NGC1961 and NGC6753, respectively. From the virial 
masses we compute the virial radii (r2oo) of the galaxies 

^ http://leda.univ-lyonl.fr 



using: 

47r 

M200 = 200pcrit yr^oo , (1) 

where Pcrit is the critical density of the Universe. We thus 
estimate that the virial radii of NGC1961 and NGC6753 
are r2oo ^ 470 kpc and r2oo ~ 440 kpc, respectively. 

3. DATA REDUCTION 
3.1. X-ray observations 

NGC1961 was observed with the European Photon 
Imaging Camera (EPIC) aboard XMM-Newton in two 
observations (ObsID: 0673170101, 0673170301) for a to- 
tal of 73.7 ks, whereas NGC6753 was observed in one 
pointing (Obsid: 0673170201) for 73.9 ks. Further de- 
tails about the observations are listed in Table [2] 

Since the main steps of th e data reduction are iden- 
tical with those outl i ned in IBourdin et al.l ()2004D and 
IBourdin fc Mazzottal (|2008D . here we only outline the 
major points. As a first step, we identify and clean 
the flare contaminated time periods for each observa- 
tion and ca mera separately, e mploying a two-step filter- 
ing process (|Nevalainen et al.. ,2005.) . The clean expo- 
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Fig. 3. — Wavelet "denoised" surface brightness image of NGC1961 based on the combined XMM-Newton EPIC datasets. The side length 
of the image is 10', which corresponds to 162.3 kpc at the distance of the galaxy. The image reveals large-scale extended emission surrounding 
the galaxy originating from the hot X-ray gas. The distribution of the hot gas is fairly symmetric. To illustrate the extent and distribution 
of the stellar light, the K-band intensity levels are overplotted. The contour levels are logarithmically scaled from 2 X lO'^ Lk,0 arcsec"^ 
to 5 X 10^ Lk arcsec"^. Within the outermost contour level, the enclosed light is ~3.2 X lO'^^ Lk q, which is ~76% of the total K-band 
lumiosity of NGC1961. 



sure time is ^30 — 80% of the original exposures. To in- 
crease the signal-to-noise ratio, we combine the data from 
the EPIC MOS and EPIC PN instruments, and we also 
merge the two pointings of NGC1961. As the goal of our 
study is to characterize the diffuse emission, bright point 
sources must be identified and removed. Therefore, we 
use the MOS data and appl y the algorithm discussed by 
IBourdin fc Mazzottal (|2008[ ) . The detected point sources 
are removed from further analysis of the diffuse emission. 

Besides the emission associated with the target galax- 
ies, a number of background components also are present. 
To account for these components, we modelled the con- 
tribution of each. The instrumental background origi- 
nates from the interactions of cosmic rays with the detec- 
tor material. The intrumental background component is 
modelled as the sum of a dual power law continuum, com- 
plemented with de tector fluorescence lin es in the 1.4 — 2 
keV energy range ()Bourdin et al.ll200l . The sky back- 
ground is modelled with three different components. The 
soft Galactic component, caused by the Local Hot Bub- 
ble, is described with two optically-thin thermal plasma 
emission models, whose temp eratures are kTi = 0.07 4 
keV and fcTa = 0.204 keV ()Kuntz fc SnowdenI [2000l ). 



The third component is a power law model with slope 
of F = 1.42, which describes the hard emission originat- 
ing from unresolv ed cosmic X-ray background sources 
(|Lumb et al.|[200l . 

Since accurate background subtraction is crucial in 
our study, we checked the accuracy of the above dis- 
cussed method by two means. First, we accounted 
for the ba ckground component s using the EPIC "blank 
sky" files (|Ca.rter fc Readll2007| ). To match the observed 
count rates, we re-normalized the background files us- 
ing the 11 — 12 keV band count rates, where the ef- 
fective area of XMM-Newton is negligible. Second, we 
analyzed the publicly available Chandra observations of 
NGC1961 (Obsid: 10528-10533). To reduce the data and 
subtract the backgr ound components, we follo wed the 
method described in iBogdan fc Gil fanovl (|2008i). Using 
these two different methods we constructed X-ray sur- 
face brightness distributions and extracted X-ray energy 
spectra of the studied galaxies (only NGC1961 for Chan- 
dra). These were compared with the profiles and spec- 
tra obtained from our standard background subtraction 
method, which models all background components. All 
three methods produce surface brightness profiles and 
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Fig. 4. — Wavelet "denoised" surface brightness image of NGC6753 based on the combined XMM-Newton EPIC datasets. The side length 
of the image is 10', which corresponds to 126.7 kpc at the distance of the galaxy. The diffuse X-ray emission, whose spectral signature is 
that of hot gas, extends beyond the optical radius of the galaxy. The overplotted K-band intensity levels demonstrate the distribution and 
extent of the stellar light. The contour levels are logarithmically scaled from 2 X lO'^ Lk,0 arcsec"^ to 5 X 10* Lk,q arcsec"^. Within the 
outermost contour level the enclosed light is ~2.8 x 10^^ Lk 0, which is ~88% of the total K-band lumiosity of the galaxy. 



X-ray spectra, which are in excellent agreement with 
each other, thereby indicating the robustness of our back- 
ground subtraction procedure. 

3.2. Near- infrared observations 

To trace the stellar light of NGC1961 and NGC6753, 
we rely on the K-band data of the Tw o-Micron All 
Sky S urvey (2MASS) Large Galaxy Atlas (jJarrett et al.l 
I2003f ). The observed background subtracted K-band 
counts (S) are converted to physical units using 

mK = KMAGZP - 2.5 log S , (2) 

where tok is the apparent K-band magnitude, and 
KMAGZP is the zero point magnitude given in the im- 
age header. The derived tuk is converted to an absolute 
magnitude and luminosity assuming that the absolute K- 
band magnitude of the Sun is Mk © = 3.28 mag. The 
total K-band luminosities of NGC1961 and NGC6753 are 
5.4 X 10^^ Lk,0 and 3.9 x 10^^ Lk,©, respectively. 

We compute the total stellar mass of NGC1961 and 
NGC6753 from their K-band luminosities and from the 
K-band mass-to- light ratios (Table [T|). The mass-to-light 
ratios are derived from the B — V color i ndices and re- 
sults of galaxy evolution modeling (|Bell et al..,2003 ). The 



mass-to-light ratios are 0.78 Mq/Lk,© for NGC1961 and 
0.81 Mq/Lk,© for NGC6753, hence their total stellar 
masses are 4.2 x 10" M© and 3.2 x 10" M©, respec- 
tively. 

3.3. Far-infrared observations 

We compute the star-formation rates (SFRs) of 
NGC1961 and NGC6753 based on their far-infrared 
emission observed by the Infrared Astronomical Satel- 
lite (IRAS), ari d applying the relation described by 
iKennicuttl (l998l ): 

SFRiR (M© yr-i) = 4.5 x IO'^^l^^, (^^g ^-i) ^ (3) 

where Lir is the total infrared luminosity in the 8 — 
1000 fim regime. To obtain t he total infrar e d lurn inosity, 
we follow the procedure of iMeurer et al.l ()1999l ) . The 
total far-infrared f lux (Ffir) in the 40 — 120 fj,m range 
(jHelou et al.lll988D is 

Ffir = 1.26 X 10-"(2.58F6o^„, -f Fioo^m) , (4) 

where Fgo/im and Fmoiim are the IRAS fluxes at 60 /zm 
and 100 fim (' Sanders et 31112003( 1. To convert the flux 
to a luminosity, we use Lfir = 47rr'|,FFiR, where D'j^ 
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Fig. 5. — Surface brightness distribution of the 0.3 — 2 keV band diffuse X-ray emission for NGC1961 (left panel) and NGC6753 (right 
panel) obtained from the combined data of XMM-Newton EPIC PN and EPIC MOS cameras. To construct the profiles we used circular 
annuli centered on the optical controid of each galaxy. Vignetting correction is applied and all background components are subtracted. 
The X-ray light distribution is compared with the K-band light profiles. The K-band profiles are normalized to match the level of X-ray 
emission in the innermost bins. 
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Fig. 6. — Surface brightness distribution of the 0.3 — 2 keV band diffuse X-ray emission for NGC1961 (left panel) and NGC6753 (right 
panel) obtained from the combined data of XMM-Newton EPIC PN and EPIC MOS cameras. The profiles are obtained from circular 
wedges, whose inner and outer radii arc 1.5' — 3.5' for NGC1961 and 1.25' — 5' for NGC6753. Vignetting correction is applied and the 
background components are subtracted. As shown in the upper x-axis, 90° corresponds to north and 180° corresponds to east. 



is the luminosit y distance of the galaxy. Since the 
iKennicutfl ()1998[ ) relation refers to the total infrared lu- 
m inosity, we conyert -L ftr to Lm applying the formula 
of iCalzetti erall (poM ): 

LiR ~ 1-75Lfir • (5) 

Although this conversion refers to the 1 — 1000 /im band 
and not to the 8 — 1000 /im band, the contribution from 
the 1 — 8 ii m range is not expec ted to exceed the few per 
cent level ([Calzetti et al.|[2000l ). yielding slightly overes- 
timated SFRs. We thus obtain SFRs of 15.5 Mq yr'^ 
in NGC1961 and 11.8 Mq yr"! in NGC6753, implying 
moderate ongoing star-formation. 

4. RESULTS 

To detect and characterize the physical properties of 
the hot X-ray coronae around NGC1961 and NGC6753, 



we employ three different techniques. First, we produce 
surface brightness images in the 0.3 — 2 keV band to 
map the spatial distribution of the diffuse X-ray emis- 
sion around the galaxies. Second, we build X-ray surface 
brightness profiles, which we compare with the stellar 
light distrubiton. Third, to measure the physical proper- 
tics of the coronal X-ray emisson, we extract their X-ray 
energy spectra. 

4.1. Images 

To map the large-scale spatial distribution of the dif- 
fuse emission around NGC1961 and NGC6753, we con- 
struct surface brightness images. The images are pro- 
duced in the 0.3 — 2 keV energy band. The particular 
choice of the energy range is advantageous because 1) the 
hot X-ray gas is predicted to have sub-keV temperatures 
(Section lO)) . hence the bulk of its emission is produced 
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Fig. 7.— X-ray energy spectra of NGC1961 (left panel) and NGC6753 (right panel) obtained from the XMM-Newton EPIC PN (top 
panel) MOSl (middle) and MOS2 (bottom) cameras. The depicted spectra are extracted for the (0.05 — 0.15)r200 radial range. The 
contribution of bright point sources are removed. The spectra are fit with an absorbed thermal plasma emission model. The black lines 
show the data set and the red lines represent the overall fit. Dark blue shows the overall background components, whereas light blue is the 
emission associated with the thermal spectrum. 



in the soft band; 2) the effective area of XMM-Newton 
detectors is the highest; and 3) the level of instrumental 
and cosmic X-ray background is fairly low in the 0.3 — 2 
keV band. 

In Figures [1] and m we show the unprocessed 0.3 — 2 keV 
band EPIC MOS images of the central 5' x 5' regions of 
NGC1961 and NGC6753, respectively. The X-ray images 
are compared with optical images from the Digitized Sky 
Survey (DSS). An apparent X-ray glow is associated with 
both galaxies within their optical extent. As discussed in 
Section [l] this X-ray emission originates from the popu- 
lation of resolved and unresolved discrete sources, addi- 
tionally truly diffuse gaseous emission also contributes. 
However, in the unprocessed images the diffuse X-ray 
emission cannot be identified beyond the optical extent 
of the galaxies, due to the faint nature of the hot X-ray 
coronae and the relatively high background level. 

In order to better visualize the faint X-ray coronae 
of the sample galaxies beyond their optical extent, we 
produce multiscale variance stabilization transform B- 
spline wavelet "dcnoiscd" surface brightness maps of the 
photon images (Starck ct al. 2009; Bourdin ct al. 2011). 
The images are background subtracted and vignetting 
corrected. To facilitate the comparison with the stellar 
light distribution, logarithmically-scaled contour levels of 
the K-band light are overplotted. 

The "denoised" 0.3 — 2 kcV band X-ray images suggest 
the presence of extended diffuse emission both around 
NGC1961 (Figure ^ and NGC6753 (Figure H]), whose 
extent significantly exceeds that of the stellar light. 
Whereas the optical light is confined to within ~1', the 
diffuse X-ray emission appears to extend out to several 
arcminutes in both galaxies. As discussed in Section [4.31 
the spectral signature of this emission is that of hot gas. 



The morphology of the X-ray emission is fairly symmet- 
ric, thereby hinting that the hot gas resides in hydro- 
static equilibrium in the dark matter halos of NGC1961 
and NGC6753. The characteristics of these X-ray coro- 
nae are further discussed in Sections 14.21 and l43l 

4.2. Profiles 

To further study the spatial distribution and the mor- 
phology of the diffuse X-ray emission around NGC1961 
and NGC6753, we construct X-ray surface brightness 
profiles. The profiles are extracted from the 0.3 — 2 keV 
energy band using the combined XMM-Newton EPIC 
data. From the profiles all background components are 
subtracted, vignetting correction is applied, and bright 
resolved point sources are masked out. 

To extract the X-ray surface brightness profiles we em- 
ploy circular annuli centered on the optical center of each 
galaxy. The distribution of the X-ray light is compared 
with the stellar light distribution, which is measured with 
the 2MASS K-band images. The K-band light profiles 
are extracted using the same circular annuli, moreover 
the locations of X-ray sources are masked out from the 
K-band images along with the population of bright fore- 
ground stars. Note that the uncertaintes in the X-ray 
surface brightness are of a statistical nature, but at large 
offsets the systematic uncertainties dominate. The pro- 
files, depicted in Figure [Sj confirm that the diffuse X- 
ray emission extends well beyond the optical radii for 
NGC1961 and NGC6753. Whereas the stehar light is 
confined to within 20 — 25 kpc (~1.5'), the diffuse X-ray 
emission is detected out to ~57 kpc (3.5') in NGC1961 
and -63 kpc (5') in NGC6753. This implies that the dif- 
fuse X-ray emission at large radii cannot be associated 
with the stellar body of the host galaxies. 
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In the wavelet "denoised" images the diffuse emission 
appears to be fairly symmetric. To further probe the 
spatial distribution of the X-ray emission, we produce 
X-ray surface brightness profiles using circular wedges 
centered on the centroid of each galaxy. Since we aim 
to probe the surface brightness beyond the optical ex- 
tent, the wedges are extracted from annuli with radii of 
1.5' - 3.5' (24.3 - 56.8 kpc) for NGC1961 and 1.25' - 5' 
(15.8 - 63.3 kpc) for NGC6753. Since the stehar light is 
negligible in this region, we do not plot the level of renor- 
malized K-band light. As before, the uncertainties are of 
a purely statistical nature. The circular profiles, shown 
in Figure HI demonstrate that the hot X-ray gas has a 
fairly uniform distribution around both galaxies at every 
position angle. This spatial distribution is morpholog- 
ically consistent with warm gas residing in hydrostatic 
equilibrium in the potential well of the galaxies, and is 
morphologically inconsistent with a startburst or AGN 
driven bipolar outflow. 

4.3. Spectra 

4.3.1. Observed parameters in the (0.05 — 0.15)r2oo region 

To study the nature and physical properties of the dif- 
fuse X-ray emission around NGC1961 and NGC6753, we 
extract X-ray energy spectra. As we aim to probe the 
hot gas beyond the optical radii, we exclude the cen- 
tral regions of the galaxies. To facilitate the compari- 
son with structure formation simulations, we convert all 
radial distances to r2oo units (Section 12. 2|) . For both 
galaxies, we extract spectra from circular annuli with 
(0.05 — 0.15)r2oo, which correspond to the physical dis- 
tances of 23.5 - 70.5 kpc for NGC1961 and 22.0 - 66.0 
kpc for NGC6753. 

In Figure [7] we show the EPIC PN and EPIC MOS 
spectra for NGC1961 and NGC6753, which demonstrate 
the presence of a soft X-ray emitting component at en- 
ergies below 2 keV. To describe this component, we em- 
ploy an absorbed thermal plasma model (WABS*APEC in 
XSPEC). When fitting the sp ectra, the column density 
is fixed at the Galactic value (|Dickey fc Lockmanlll990l ) 
and the abundance is left fre e to vary. For the spectra l 
fitting the abundance table of iGrevesse fc Sauvall ()1998f ) 
is used. The thermal components describe the spectra 
well, demonstrating that the diffuse X-ray emission orig- 
inates from hot gas. 

Based on the best-fit spectra we measure the average 
parameters of the hot gas within the (0.05 — 0.15)r2oo 
region. For both NGC1961 and NGC6753, we obtained 
fairly similar gas temperatures {kT ^ 0.6 — 0.7 keV) and 
sub-solar abundances (0.12 — 0.13 Solar). From the emis- 
sion measure of the spectra (J nenndV), we compute the 
average electron number density, the confined gas mass, 
and the average cooling time. As a caveat we mention 
that the emission measure and abundances of the thermal 
component are strongly anticorrelated (Section 16. 3p . To 
compute the hot gas mass, we assumed a spherically sym- 
mnictric gas distribution, constant density, and a volume 
filling factor of unity. Given the volume of the studied 
regions, we obtained hot gas masses of about ~10^° Mq 
and electron number densities of (3.5 — 4) x 10"'* cm~'^ for 
both galaxies. The cooling time of the gas is computed 
from tcooi = (3fcr)/(neA(T)), where A(T) is the cool- 
ing function at the best-fit temperature and abundance. 



For both NGC1961 and NGC6753 we obtain tcooi 50 
Gyrs, which is many times the Hubble-time. The un- 
certainties of the emission measure (and other quantities 
derived from its value) are estimated from the parameter 
covariance matrix by assuming a quadratic distribu- 
tion. The parameters of the warm gas are listed in Table 

m 

From the best-fit spectra, we also compute the X-ray 
luminosity of the hot gas in the (0.05 — 0.15)r2oo radial 
cut. In the 0.5 — 2 keV band we obtain the absorbed lu- 
minosities of io.5-2kov,abs — (2.0 ±0.6) X lO'*^ erg s~* for 
NGC1961 and io.5-2kcV,abs = (2.5 ± 0.4) x 10^° erg s-^ 
Using the Galactic column densities (|Dickev fc LockmanI 
119901 ) we convert these values to unabsorbed luminosi- 
ties. Finally, we also derive the corresponding bolo- 
metric luminosities of the X-ray coronae, and we find 
Lboi = (5.8 ±1.7) x 10-*" erg s^* and Lboi = (6.3 ±0.9) x 
10-*° erg for NGC1961 and NGC6753, respectively. 
The obtained fluxes and luminosities are listed in Table 

m 

4.3.2. Upper limits in the (0.15 — 0.30)r2oo region 

In the absence of hot gas detection beyond ~0. 15^200, 
we place upper limits on its physical parameters in the 
(0.15 — 0.30)r2oo region. This region is well within 
the XMM-Newton field-of-view, and corresponds to the 
radial range of 70.5 - 141.0 kpc for NGC1961 and 
66.0 — 132.0 kpc for NGC6753. From the observed count 
rates within these regions we derive 3(T upper limits us- 
ing the EPIC PN data. To constrain the gas proper- 
ties, we rely on the gas parameters observed within the 
(0.05 — 0.15)r2oo region. In particular, we assume the 
abundance of 0.12 Solar for NGC1961 and 0.13 Solar for 
NGC6753, and we adopt kT — 0.6 keV for both galaxies. 

Using the appropriate response files, we derive up- 
per limits on the emission measure, which we use to 
place upper limits on the confined gas mass, the elec- 
tron number density, and the cooling time. Within the 
(0.15 — 0.30)r2oo region, the 3cr upper limit on the gas 
mass is < 1.6 x 10*" Mq and < 1.3 x 10*° Mq for 
NGC1961 and NGC6753, respectively. For the average 
electron number density we derive < 6.4 x 10~^ cm~'^ 
in NGC1961 and < 6.2 x 10"^ cm'^ in NGC6753. The 
cooling times for both galaxies are longer than 250 Gyrs, 
hence significantly exceed the Hubble-time, implying the 
quasi-static nature of the hot gas. The obtained upper 
limits are listed in Tabled 

Based on the upper limits on the emission measure, 
we also constrain the X-ray luminosities within the 
(0.15 - 0.30)r2oo region (Table S]). The upper hmits 
on the absorbed 0.5 — 2 keV band luminosities are < 
4.8 X 10^9 erg s"* and < 2.5 x lO^^ erg s'^ for NGC1961 
and NGC6753, respectively. These values correspond to 
the bolometric luminosities of < 1.4 x 10'*° erg s^* in 
NGC1961 and < 6.5 x lO^^ erg s"* in NGC6753. 

5. COMPARISON WITH NUMERICAL SIMULATIONS 

5.1. Numerical codes 

In this paper the observed properties of X-ray coro- 
nae around NGC1961 and NGC6753 are confronted with 
hydrodynamical simulations of structure formation. In 
particular, we probe the recently developed moving mesh 
code AREPO (|Springel et al., ,201ft ) and the traditional 
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TABLE 3 

The parameters of the extended hot gaseous coronae around NGC1961 and NGC6753. 



Galaxy 


Radial range 




kT 


Abundance 


Mgas 


ric 


^cool 




''200 


_ ■a 

cm 


keV 




Me 


cm ^ 


Gyr 


NGC1961 


0.05 - 0.15 


(5.2 ± 1.5) X 10'33 




0.12 ±0.03 


(1.2 ±0.2) X 10^° 


(3.5 ±0.5) X 10"* 


45 ±6 


NGC6753 


0.05 - 0.15 


(5.3 ±0.8) X 1063 


69+" '"^ 
'-'■"^-0.07 


0.13 ±0.02 


(1.1 ± 0.1) X 10^° 


(3.9 ±0.3) X 10-4 


51 ±4 


NGC1961 


0.15 - 0.30 


< 1.2 X 10*^3 


0.6t 


0.121' 


< 1.6 X lO^o 


< 6.4 X 10-5 


> 258 


NGC6753 


0.15 - 0.30 


< 9.6 X lO'*^ 


0.6t 


0.13t 


< 1.3 X lO^o 


< 6.2 X 10-5 


> 266 



1 These parameters were fixed at the given value. 



TABLE 4 

The flux and luminosity of the extended hot gaseous coronae around NGC1961 and NGC6753. 



Galaxy 


Radial range 


-Fb.5-2kcV,abs 


-f'0.5-2kcV,abs 


-^'0.5-2kcV,unabs 


Lhoi 




''200 


erg cm" 


erg s"^ 


erg 


erg 


NGC1961 


0.05 - 0.15 


(5.4 ±1.6) X 10-" 


(2.0 ±0.6) X lO-*" 


(2.8 ±0.8) X 10'"' 


(5.8 ± 1.7) X 10"° 


NGC6753 


0.05 - 0.15 


(1.1 ±0.2) X 10"" 


(2.5 ±0.4) X 10''° 


(3.1 ±0.5) X 10"" 


(6.3 ±0.9) X 10"" 


NGC1961 


0.15 - 0.30 


< 1.3 X 10-" 


< 4.8 X 10^9 


< 6.7 X lO^'' 


< 1.4 X 10"" 


NGC6753 


0.15 - 0.30 


< 1.1 X 10-" 


< 2.5 X lO^^* 


< 3.1 X 10^9 


< 6.5 X 1039 



SPH-based code GADGET ("Springcll I2005D . The codes 
and the methodol ogy of the simulations are discussed in 
full particulars in iVogelsberger et al.l ()2012D . 

GADGET is a widely used and well tested SPH code, 
which simultaneously conserves energy and entropy 
despite the use of fully adaptive smoothing lengths 
(jSpringel fc Hernou ist 2002). The gravity calculation is 
split into short-range and long-range components, where 
the short-range forces are calculated with a hierarchi- 
cal q ct-tree algorithm (jBarnes fc HutI [l986l : IHernauistl 
I1987I) and the long-range forces are evaluated with a 
parti cle mesh (PM) method (e.g. IHocknev fc EastwoodI 
|1981|) . The simulations presented here are based on a 
default setting of the n umerical SPH parameters (see 
IVogelsberger et al]|2012l for details) . AREPO is a second- 
order accurate finite volume method that solves the Eu- 
ler equations using piece-wise linear reconstruction and 
a calculation of hydrodynamical fluxes at each cell inter- 
face with an exact Riemann solver. The basic solution 
strategy of the code is that of the well-known MUSCL- 
Hancock scheme. AREPO employs an unstructured mesh 
based on a Voronoi tessellation using a set of mesh- 
generating points. This mesh is allowed to move freely as 
a function of time. In the simulations presented here, the 
motion of the mesh-generating points is tied to the hy- 
drodynamical flow resulting in a quasi-Lagrangian char- 
acter and makes the mesh automatically adaptive. Such 
a scheme combines the advantages of classical adaptive 
mesh refinement (AMR) and SPH codes. For example, 
shocks, contact discontinuities, and mixing are signifi- 
cantly better resolved in AREPO than in classical SPH 
codes like gadget. 

The simulations discussed here adopt the cosmological 
parameters fim.o = 0.27, f^A^o = 0.73, f2b,o = 0.045, 
as = 0.8, ris = 0.95 and Hq = 100 /ikms^i Mpc"^ = 
70kms~^ Mpc~^. These parameters are consistent with 
the most recent WMAP-7 measurements (jKomatsu et al.l 
1201 If) . The simulation covers a volume of (20 Mpc)'^ 
sampled by 2 x 512^^ dark matter particles/cells with a 
gravitational softening length of 1 h^^ kpc comoving 



The gadget and AREPO simulations follow the same 
physics. In particular, the simulations consist of a colli- 
sionless dark matter fluid and an ideal baryonic gas aug- 
mented with additional terms that account for radiative 
processes (optically-thin radiative cooling of a primor- 
dial mixture of helium and hydrogen plus a uniform, 
time-dependent ionising UV background) and star for- 
mation. Gas can collapse to high density and turn into 
stars. Both codes describe this process with the star 
formation and supernova fe edback model introduced in 
iSpringel fc HernquistI (|2003l ) assuming a Salpeter IMF. 

A comparison between the moving mesh code (arepo) 
and the SPH code (gadget) has been perforrned re- 
cently in a series of papers (iVogelsberger et al.l 



Keres et al 



1^511 ISiiacki et al.N2012l: 



Bird et all (20121) . In their study. 



Torrev et al. 



2012 



2011 



Vogelsberger et al] 



(|2012f) recognized a number of important differences be 
tween the two codes. Most importantly, significantly 
more gas cools out of halos at low redshifts in arepo than 
in GADGET, thereby yielding more massive and more ex- 
tended stellar disks. This difference can be attributed 
to the stronger cooling flows that evolve in arepo at 
low redshifts, yielding larger star-formation rates. At 
the same time, however, the hot X-ray coronae around 
massive galaxies are less dense, hence less luminous and 
have lower mass in arepo than in gadget. Addition- 
ally, the luminosity weighted gas temperatures are some- 
what lower in AREPO than in gadget. The hotter gas 
temperatures and the weaker gas cooling in the GAD- 
GET code can be explained by the different dissipative 
heating rates. Indeed, arepo produces higher dissipa- 
tive heating in the halo infall regions, and gadget pro- 
duces higher dissip ation in the outer regions of halos. 
IVogelsberger et all (2012) argue that the higher dissipa- 
tion in_GADGET_Js_ likely to be of spurious nature (see 
also iBauer fc Spring cl 2012). Thus, the differences de- 
scribed above between the arepo and gadget simula- 
tions are purely due to the inaccuracies of the hydrody- 
namic solver. 

Given the above differences between the arepo and 
GADGET simulations, it is essential to observationally test 
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Fig. 8. — Predicted bolometric X-ray luminosities of the X-ray coronae around disk galaxies as a function of stellar mass for the 
(0.05 — 0.15)r'2oo (left panel) and (0.15 — 0.30)r2oo radial cuts (right panel). Galaxies simulated by AREPO are shown with small boxes, 
whereas galaxies simulated by GADGET are marked with small circles. The observed X-ray luminosities and upper limits for NGC1961 and 
NGC6753 are shown with the big symbols. The statistical uncertainties on the observed X-ray luminosities arc comparable with the size 
of the symbols. 
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Fig. 9. — Predicted hot gas mass of the X-ray coronae around disk galaxies as a function of stellar mass for the (0.05 — 0.15)r200 (left 
panel) and (0.15 — 0.30)r2oo regions (right panel). Only gas particles with temperature above 0.1 keV were used to compute the gas mass 
in the simulations. Galaxies simulated by the AREPO code are shown with small boxes, whereas galaxies simulated by the GADGET code are 
marked with small circles. The measured hot gas masses and upper limits for NGC1961 and NGC6753 are shown with the big symbols. 
The statistical uncertainties on the hot gas masses are comparable with the size of the symbols. 



these models. Confronting the observed and predicted 
properties of X-ray coronae around massive spiral galax- 
ies is advantageous for the following two reasons. First, 
the existence of hot gaseous coronae around massive disk 
gala xies is a ba sic prediction of galaxy formation mod- 
els (jWhite fc R ecs 1978), hence disk galaxies, such as 
NGC1961 and NGC6753, offer a clean observational test 
of these models. Second, the extended X-ray coronae 
are well-resolved in numerical simulations, hence their 
evolution can be followed through a large redshift range. 
Thus, in the following sections we compare the observed 
properties of hot gasesous coronae around NGC1961 and 
NGC6753 with those predicted by the AREPO and the 
GADGET simulations. 

5.2. X-ray luminosities 



As a first probe, we confront the observed X-ray lu- 
minosities around NGC1961 and NGC6753 with those 
predicted by the AREPO and GADGET simulations. To fa- 
cilitate the comparison with present-day and future sim- 
ulations, all radial distances are converted to r2oo units, 
where r2oo is defined as the virial radius. The r2oo radii 
of NGC1961 and NGC6753 are -470 kpc and -440 kpc, 
respectively (Section I2.2[) . We note that the simulated 
galaxies are not members of rich galaxy clusters. 

For both NGC1961 and NGC6753, we compare the 
observed and predicted bolometric X-ray luminosities in 
two circular annuli with radii of (0.05 — 0.15)r2oo (inner 
region) and (0.15 — 0.30)r2oo (outer region). Note that 
the central regions (< 0.057-200) are excluded due to the 
various contaminating factors associated with the stellar 
body of the galaxies (SectionH]). In the (0.05 — 0.15)r2oo 
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Fig. 10. — Observed and predicted electron density profiles as the function of the radius for the GADGET (left panel) and AREPO (right 
panel) simulations. Note that the radius is measured in units of r2oo- The thick solid lines show the obs erve d density profiles, whereas the 
short dashed thick lines show the extrapolated density profiles based on the best-fit /3-models (Section 15.31 1. The statistical uncertainties 
on the observed profiles are comparable with the difference between the two profiles. The predicted density profiles originate from five 
representative simulated galaxies, whose stellar masses are comparable with those of NGC1961 and NGC6753. 



TABLE 5 

The best fit parameters of hot gas density profiles. 



Galaxy 


no 




rl 


a 


/3 


ft 


7t 




(cm"-^) 


(kpc) 


(kpc) 










NGC1961 


1.0 X 10"^ 


18.1 


400 


4.40 


0.29 





3.0 


NGC6753 


1.8 X 10"^ 


9.7 


400 


6.04 


0.36 





3.0 



t These parameters were fixed at the given value. 

region, the X-ray luminosities are obtained from the best- 
fit spectra. In the (0.15 — 0.30)r2oo region, in the absence 
of X-ray gas detection, we place 3cr upper limits on the 
luminosities using adopted gas parameters (Section 14. 3p . 
The bolometric X-ray luminosities and upper limits for 
NGC1961 and NGC6753 are listed in Table H 

In Figure |8] we show the predicted bolometric X-ray 
luminosites as a function of the stellar mass for the 
(0.05 — 0.15)r2oo and (0.15 — 0.30)r2oo regions. In these 
plots we mark the corresponding observed values and 
upper limits for NGC1961 and NGC6753. It is appar- 
ent from Figure [8] that, on average, the predicted X-ray 
coronae are more luminous in GADGET than in AREPO. 
Although this difference is clear in the inner region, it 
becomes even more evident in the outer regions. As a 
caveat, we note that above ~2 x 10^^ Mq this difference 
is less evident, especially in the (0.05 — 0.15)r2oo radial 
cut. In the inner region, the observed luminosities for 
NGC1961 and NGC6753 fall closer to the gadget pre- 
dictions. However, given the relatively large scatter at 
the high mass end, no strong conclusions can be drawn. 
In the outer region, the 3cr upper limits are approxi- 
mately consistent with the AREPO model predictions, but 
fall a factor of ^2 — 10 short of the gadget predictions. 

5.3. Hot gas mass and density profiles 

In Figure [S] we present the predicted gas mass in the 
(0.05 — 0.15)r2oo and (0.15 — 0.30)r2oo regions as a func- 
tion of the stellar mass. In both AREPO and GADGET 
simulations we only consider the hot gas components, 
whose temperature is above 0.1 keV. The observed gas 



masses and upper limits are given in Table [31 

Figure [S] demonstrates that gadget predicts signifi- 
cantly more massive X-ray coronae than AREPO. This 
difference is significantly larger in the outer region, 
thereby hinting that the density profiles of the hot gas 
coronae fall more rapidly in the AREPO than in GADGET 
simulations. In the inner region, the observed gas mass 
is somewhat higher than that predicted by the moving 
mesh code but falls short of the SPH simulations. In the 
outer region, the 3a upper limits on the gas mass are 
consistent with the AREPO predictions but are factors of 
2 — 10 less than expected from the GADGET simulations. 

To further probe the structure formation models, we 
compare the observed and predicted electron density pro- 
files. Since the hot gas is only detected out to '^0.15r2oo, 
we extrapolate our measurements out to the virial radius, 
which enables a more detailed comparison with the sim- 
ulations. To this end, we describe the observed surface 
brightness distribution with a modified /3-model. Under 
the assumption of constant metallicity gas, we derive the 
em ission measure p r ofile a t a projected radius (r) follow- 
ing |VikhiininitaLl ()2006[ ): 



-p-e - -0 ^ ^2/^2)3;9-a/2 + ^7 /r]Yh ' ' 

This /3-model consists of a power law-type cusp and also 
describes the possible steepening of X-ray surface bright- 
ness profiles at large radii ( Vikh linin et a l. 2006). The 
obtained best-fit parameters of the emission measure pro- 
file are listed in Table \E\ To compare the observed elec- 
tron density profiles with those produced by AREPO and 
GADGET, we extract the density profiles of five represen- 
tative simulated galaxies, whose stellar masses are com- 
parable to NGC1961 and NGC6753. 

The comparison between observed and simulated elec- 
tron density profiles are depicted in Figure [101 Although 
the density profiles predicted by the AREPO simulations 
exhibit a large scatter, they are in good agreement within 
~0.4r2oo (left panel of Figure [TO\i with the profiles of 
NGC1961 and NGC6753. Indeed, not only the shape 
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FlG. 11. — Predicted hot gas temperatures of the X-ray coronae 
as a function of the stellar mass for the (0.05 — 0.15)r200 range. 
Simulated galaxies are shown with the small boxes (arepo) and 
small circles (gadget). The observed hot gas temperatures for 
NGC1961 and NGC6753, obtained from the best-fit spectra, are 
shown with the big filled symbols. 
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Fig. 12. — Predicted abundances of the X-ray coronae vs stel- 
lar mass for the (0.05 — 0.15)r200 range. Galaxies simulated by 
AREPO are shown with boxes, and galaxies simulated by GADGET 
are marked with circles. The observed abundances for NGC1961 
and NGC6753 are shown with the big symbols. The statistical un- 
certainties on the measured values are comparable with the size of 
the symbols. 



of the profiles, but also the observed and predicted rie 
values are comparable. We note, however, that within 
'-~^0.15r2oo the observed densities are at the high end of 
the model predictions, hence in the (0.05 — 0.15)r2oo re- 
gion the observed X-ray gas masses (and luminosities) are 
somewhat higher than those expected by the AREPO sim- 
ulations. As a caveat, we mention that within ~0.05r2oo 
the X-ray emission is dominated by gas produced in the 
course of the stellar evolution of the galaxies, therefore 
this region does not test the extended X-ray coronae. Be- 
yond ~0.4r2oo the density drops sharper than the extrap- 
olated profiles of NGC1961 and NGC6753. The trun- 
cated density profile in the simulation is most likely due 
to the absence of efficient supernova and AGN feedback 
in the models. Implementing efficient energy feedback 
will likely result in more extended gas distributions, re- 
sulting in a better agreement with the observed profiles 
at large galactocentric radii. 

The density profiles predicted by the gadget sim- 
ulations exhibit strikingly different distributions than 
the observed ones at all radii (right panel of Figure 
[TU|) . Indeed, the gadget profiles do not only predict 
too high densities in the (0.05 — 0.4)r2oo radial range, 
but within this region the simulated density profiles are 
significantly flatter than those observed for NGC1961 
and NGC6753. Although at large radu (> 0.4r2oo) the 
predicted and observed (extrapolated) electron densities 
agree with each other, the shape of the profiles signifi- 
cantly deviate from each other. Indeed, the simulated 
profiles drop much faster than the extrapolated ones at 
radii beyond > 0.4r2oo, which, at least partly, may be ex- 
plained by the absence of efficient supernova and AGN 
feedback in the model. 

Based on the density profiles we conclude that (1) 
AREPO describes the observed density profiles better than 
GADGET at radii < 0.4r2oo; and (2) at larger radii both 
simulations fail to reproduce the extrapolated profiles, 
presumably due to the absence of efficient feedback in 
the models. 



10' 



10' 



10' 



10' 




OO CD 



AREPO 
GADGET 



10 



10. 



10 



Fig. 13. — Virial mass as a function of stellar mass for the AREPO 
and GADGET simulations. Galaxies simulated by the AREPO code 
are shown with boxes, whereas galaxies simulated by the GADGET 
code are marked with circles. Note that both simulations suffer by 
the so-called overcooling problem, that is too much gas cools out 
of the halos, thereby producing too massive stellar disks relative 
to their dark matter halos. The stellar masses of NGC1961 and 
NGC6753 are 4.2 x lO" Mq and 3.2 X lO^^ Mq, respectively. 



5.4. Hot gas temperature 

In Figure [TT] we depict the predicted hot gas temper- 
atures in the (0.05 — 0.15)r2oo radial range as the func- 
tion of the stellar mass. The predicted temperatures are 
luminosity- weighted values. 

In general, GADGET predicts higher gas temperatures 
than AREPO throughout the entire stellar mass range. 
For both simulations, the expected temperatures are in 
the relatively narrow kT = 0.3 — 0.7 keV range for galax- 
ies with ~4 X 10" Mq. For NGC1961 and NGC6753 
the observed gas temperatures are feT ^ 0.6 — 0.7 keV, 
which places them at the high end of the predictions. 
The observed values are in fair agreement with the SPH 
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simulation and exceed the values predicted by the mov- 
ing mesh code. 

However, at present, neither simulation code imple- 
ment efficient energy feedback from supernovae and su- 
permassive black holes, which serves as an additional 
heating source. Therefore, including the efficient feed- 
back in the simulations will yield higher gas tempera- 
tures, and possibly better agreement with the observed 
values. We note that based on the presently available 
X-ray data, the temperature distribution of the hot gas 
cannot be studied in detail. To better explore the tem- 
perature structure of X-ray coronae, further deep obser- 
vations are required. 

5.5. Abundances of the hot gas 

As a final test, we investigate the observed and pre- 
dicted abundances of the hot X-ray coronae in the 
(0.05 — 0.15)r2oo radial cut. The observed abundances 
for NGC1961 and NGC6753 are the best-fit values from 
the spectral fitting procedure, which are given in Table 
[3l The predicted values are luminosity-weighted abun- 
dances. 

In Figure [12] we show the predicted abundances for the 
AREPO and GADGET simulations as a function of stellar 
mass. The observed values for NGC1961 and NGC6753 
are overplotted. For the two structure formation simu- 
lations, there is a striking difference between the pre- 
dicted abundances. Whereas for galaxies with a few 
xlO^-^ M© stellar mass AREPO predicts ■^O.l Solar abun- 
dance. The expected values in gadget are about two 
orders of magnitude less, ~10~^ Solar. The observed val- 
ues for NGC1961 and NGC6753 are ~0.1 Solar, which are 
in good agreement with that predicted by AREPO. We 
thus conclude that the observed abundances of the hot 
gas around NGC1961 and NGC6753 strongly favor the 
moving mesh code over the SPH simulations. 

The abundance difference between the codes is caused 
by two effects. First, due to the increased cooling, AREPO 
exhibits higher star formation rates than GADGET, which 
is attributed to the different dissipative heating rates 
(Section 15. 1|) . Since metal production in the simula- 
tions is tied to star formation assuming a fixed yield, 
the higher star formation in AREPO automatically leads 
to a production of higher abundances. Second, the mix- 
ing is more efficient in AREPO than in gadget. This 
point is related to the intrinsic mixing problems of stan- 
dard SPH codes. In the gadget simulations, metals 
are transported by individual SPH particles and they do 
not undergo any explicit metal diffusion process. Con- 
trarily, the metals are advected and can mix well within 
the gas in the AREPO code. Accordingly, metals stick 
more closely to the central part of galaxies and haloes in 
GADGET compared to AREPO. This difference becomes 
particularly apparent for more massive galaxies, where 
enrichment of halos through mixing is highly ineffective 
in GADGET (Figure [T^ . In massive galaxies the abun- 
dances differ by about two orders of magnitude between 
the two simulations, but this difference gets progressively 
lower in smaller mass systems. Finally, we note that the 
inclusion of explicit supernova winds is another source of 
halo gas enrichment, which is not included in our simu- 
lations. 

5.6. Limitations of the present simulations 



Although AGN feedback and galactic winds are 
believed to have a profound impact on the evolu- 
tion of galaxies, the galaxy formation simulations dis- 
cussed in this paper do not incorporate efficient en- 
ergy feedback from AGN and supernovae. The im- 
pact of AGN feedback is observed i n different scales 
in the local Universe : in galaxies ([Finoguenov et al.l 
l200ll:lJones et al.ll2002l: [Kraft et a l.ll2011'). galaxy group s 
pavid et al.l2009tlGitti et al.ll201 0: Randall e t al.l201lD 
and galaxy clusters (Fabian et al. 2000; For man et al.l 
J005; Nulscn ct al. 2005; Bogdan ct al. 201 id). More- 
over, theoretical studies also point out that the AGN 
feedback is capable to heat the cool gas and expel it 
to larger radii, thereby quenching the star-formation 
and regulating the gro wth of supermassiye black holes 
J^ilk & Rccs 1998; Di Matteo et all [20051 iSiiacki et al.1 
120071) . Galactic winds driven by the energy input of su- 
pernovae eject metal-rich gas out to large radii, thereby 
enric hing the extraplanar gas (jOppenheimer fc Pavel 
|2006[) . Additionally, galactic winds play a major role 
in reprod ucing the cosmic sta r-formation history at high 
redshifts (jSchave et al.l[20Tol ). Furthermore, implement- 
ing efficient AGN and supernova feedback will address 
one of the major shortcomings of the present models: 
the so-called overcooling of galaxies. Without efficient 
feedback too much gas cools out of the halos both in 
the AREPO and the GADGET simulations, hence produce 
overly massive stellar components relative to their dark 
matter halos. Although both simulation codes are af- 
fected by the overcooling problem, galaxies in AREPO are 
more massive than their counterparts in gadget (Figure 

m- 

The present models also do not incorporate the con- 
tribution of metal cooling to the emission. However, the 
cooling of metal species plays a significant role at the ob - 
served sub-keV temperatures (jvan de Voort et"aIll201lD . 
Indeed, metal cooling notably influences the tempera- 
ture structure, hence the X-ray luminosity of the hot gas. 
Metal line cooling also increases the star formation rates, 
and affects the build-up of galaxies. On the other hand, 
AGN and supernova feedback typically heat the gas and 
decrease the star formation rates, such that missing feed- 
back and missing metal line cooling partially cancel each 
other. 

In subsequent papers we will investigate structure for- 
mation simulations, which incorporate the above dis- 
cussed missing physics. Thus, the forthcoming studies 
will lead to a physically more realistic comparisons be- 
tween the observed X-ray coronae and numerical models. 

5.7. Synopsis about the 
AREPO and GADGET simulations 

In this section we compared the observed physical 
properties of hot X-ray coronae around two massive spi- 
ral galaxies with those predicted by two structure forma- 
tion simulations. The moving mesh code (arepo) and 
the SPH based code (gadget) uses fully identical sub- 
resolution physics and gravity solver, but they employ 
different methods to solve the equations of hydrodynam- 
ics. Based on the comparison we obtained the following 
major conclusions: 

1. In the (0.05 — 0.15)r2oo radial range the X-ray lu- 
minosities and hot gas masses broadly agree with 
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both simulations. However, in the (0.15 — 0.30)r2oo 
region the upper limits for X-ray luminosities and 
gas masses are more consistent with AREPO, since 
GADGET overpredicts both quantities. 

2. The observed (and partly extrapolated) density 
profiles in the (0.05— 0.4)r2oo radial range are in fair 
agreement with those predicted by AREPO. Con- 
trarily, GADGET predicts markedly different density 
profiles, which decline much slower and disagree 
with the shape of the observed profiles at all radii. 

3. The observed abundances in the (0.05 — 0.15)r2oo 
region are in good agreement with those predicted 
by AREPO, but exceed the predictions of gadget 
by about two orders of magnitude. 

Although the above listed points suggest that the mov- 
ing mesh code gives a significantly better description of 
structure formation, further observational and theoret- 
ical efforts are required before definite conclusions can 
be drawn. On the observational side it is necessary to 
further push the limits of X-ray observations, and probe 
the X-ray coronae at even larger radii. Such observations 
would allow the more precise comparison of observations 
and numerical models. Moreover, the parameter space 
needs to be better explored, that is the X-ray coronae of 
other, preferably lower mass, galaxies should be observed 
and detected. On the theoretical side, efficient energy 
feedback from AGN and supernovae must be included 
along with cooling from metal species. The comparison 
with physically more realistic models will lead to a sig- 
nificantly better understanding of the physical processes, 
which influence galaxy formation models. 

6. DISCUSSION 
6.1. The role of starburst driven winds 

Although the goal of the present work is to ex- 
plore the hot gaseous X-ray coronae around massive 
spiral galaxies, it is feasible that certain fraction of 
the detected hot gas was uplifted from the galactic 
disks. In principle, galaxies with sufficiently high spe- 
cific star-formation rates are capable of driving star- 
burst dri ven winds. Based on theoretical computa- 
tions, Strickland et al.l ([2004b ) estimated that supernova 
blowout occurs at the (surface area) specific SN rate 
>25 SN Myr~^ kpc~^, where the surface area is derived 
as the square of the D25 diamet er. In fair agreem e nt with 
the theoretical computations, [Strickland et al.l (l2004bl ) 
detects extraplanar X-ray emission in galaxies, whose 
specific SN rates are >40 SN Myr'^ kpc"^. 

To probe the importance of SN driven winds in 
NGC1961 and NGC 6753, we compu t e thei r area specific 
SN rates. Following iHeckman etall ()1990[) . the SN rate 
is derived from i?gN = 0.2Lfir/10^^ Lq, where Lpm, 
is the total far-infrared luminosity. For NGC1961 and 
NGC6753 we observe iFiR,NGCi96i = 9.0 x 10^*^ L© and 
■C'FiR,NGC6753 = 6.8 X 10^° Lq fSection lS^ . which values 
correspond to the SN rates of i?SN,NGCi96i = 0.18 yr~^ 
and i?sN,NGC6753 = 0.14 yr~^, respectively. Based on 
these values and the D25 major axis diameters, we es- 
timate the area specific SN rates of ~36 Myr~^ kpc~^ 
for NGC1961 and ~97 Myr^^ kpc'^ for NGC6753. The 



value obtained for NGC1961 is comparable albeit some- 
what lower than the observational lower limit required 
for SN blowout. In NGC6753 the specific SN rate is 
higher than the critical rate for SN blowout, hinting that 
fraction of the detected warm gas may be uplfited from 
the disk. 

Despite the relatively high specific SN rates, the fol- 
lowing four arguments suggest that the detected X-ray 
coronae around NGC1961 and NGC6753 cannot be en- 
tirely attributed to starburst driven winds. 1) The mor- 
phology of the X-ray gas (Figures [3l IH [6]) is fairly sym- 
metric. This indicates that the hot gas resides in hy- 
drostatic equilibrium in the potential well of the galax- 
ies, and is not in a bipolar starburst driven outflow. 
2) Although the specific SN rates are comparable or 
exceed the criteria for SN blowout, the derived values 
are signifcantly lower than that observed in starburst 
galaxies, for example in M82 (~1200 M q Myr~^ kpc^^) 
or NG C253 (~180 Mq Myr"! kpc^^) (jstrickland et al.l 
(IomO). 3) The extraplanar X-ray emis sion in starburst 
galaxies is confined to within ^15 kpc (jStrickland et al.l 
|2004a), whereas we detect hot gas out to ^60 kpc around 
NGC1961 and NGC6753. Additionally, we emphasize 
that the central 0.05r2oo ('--"23 kpc) regions are excluded 
when studying the properties of the diffuse X-ray emis- 
sion. 4) The cooling time of the hot gas is on the time 
scale of ~50 Gyrs in the inner region. Therefore, the 
hot gas can be treated as quasi-static gas as opposed 
to a starburst driven wind. Taken together, these evi- 
dences suggest that only fraction of the warm gas around 
NGC1961 and NGC6753 could originate from starburst 
driven winds. 

The above considerations demonstrate that at the 
present epoch the X-ray appearance of NGC1961 and 
NGC6753 is not determined by starburts driven winds. 
However, at higher redshifts gas-rich galaxy mergers oc- 
cured more frequently, which could have triggered more 
active phases, possibly leading to starburst driven winds. 
Such winds do not only expel a notable amount of gas 
to larger radii, but also play a major role in enriching 
the hot extraplanar gas, which is reflected by the non- 
zero abundance of the hot gas around NGC1961 and 
NGC6753. Thus, at some point of their evolution, the 
hot X-ray coronae are likely to be influenced by starburst 
driven winds. 

6.2. Baryon mass fraction 

The detection of hot X-ray coronae around NGC1961 
and NGG6753 allows us to estimate their baryon budget 
and derive their baryon mass fraction. To account for 
the baryon content, we consider the total hot X-ray gas 
mass, the stellar mass of the galaxies, the cold gas mass, 
and the stellar mass of the group member galaxies. The 
detailed mass budget is given in Table |6l 

In massive galaxies, such as NGC1961 and NGC6753, 
the hot X-ray emitting gas adds a notable contribution 
to the baryon budget. Therefore, we derive the confined 
hot gas mass by extrapolating the density profiles out 
to the virial radius (Section 15. 3p . For NGC1961 and 
NGC6753 the obtained total hot gas masses within r2oo 
are 9.4 x 10" Mq for NGC1961 and 4.3 x 10" M©, re- 
spectively. We emphasize that these numbers should only 
be considered as a crude estimate due to the following 
systematical and statistical uncertainties. 1) The ma- 
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TABLE 6 

Mass budget of NGC1961 and NGC6753. 

Galaxy Af* M'gas.hot A^gas.cold Mgroup Afb.tot My)M fh 

(Mq) (Mq) (Mq) (Mq) (Mq) (Mq) 

(1) (2) (3) (4) (5) (6) (7) 

NGC1961 4.2 X lO" 9.4 x lO" 4.7 x lO^o 1.4 x 10" 1.55 x lO^^ 1.2 x lO^^ 0.11 

NGC6753 3.2 x lO" 4.3 x lO" 2.0 x lO^" 2.7 x lO" 1.04 x lO^^ 1.0 x 10^^ 0.09 



Note. Columns are as follows. (1) Total stellar mass of th e sa mple galaxies. (2) Estimated hot X-ray gas mass within the virial radius, 
obtained from the extrapolation of the /3— models (Section 16.21 1. (3) Cold gas mass. (4) Total stellar mass of group members. (5) Total 
baryon mass, which is the sum of columns (1) — (4). (6) Estimated mass of the dark matter halo. (7) Baryon mass fraction, defined as 
fh = Mb,tot/(A/DM + Afb,tot)- 



jor source of uncertainty originates from the unexplored 
nature of the hot gas beyond 0.15r2oo radius. Indeed, 
a temperature and/or abundance gradient could drasti- 
cally alter the extrapolated density profiles hence the to- 
tal gas mass within the virial radius. 2) The hot gas may 
not have a spherically symmetric distribution at large 
radii. 3) The applied r2oo radii may deviate from the 
real r2oo radii. 4) Statistical uncertainties in the best-fit 
parameters of the /3-model. 

The stellar mass is derived from the 2MASS K-band 
images and the K-band mass-to-light ratios (Section 
13. 2|) . For the cold gas mas s we use the HI mass of 
4.7 X IQi" Mq for NGC1961 (IHaan et al.ll2n08h . and we 
estimate the HI mass of about 2 x lO^" Mq for NGC6753 
nes fc GiovanellilfTosl. Since both ga laxies are lo- 
cated in poor groups (fFouque et al.l 11993 ). the baryon 
mass of these galaxies should also be included in the 
baryon budget. The baryon mass of the group members 
is dominated by their stellar mass, which we estimate 
based on their K-band luminosity, using M^,/Lk — 0.8, 
and assuming that they lie at the same distance as 
NGC1961 and NGC6753. 

The total baryonic mass of NGC1961 and NGC6753 
areMb,tot = l.SSxlO^^ jyj^ andMb,tot = 1.04x10^2 Mq, 
respectively (Table [5]). The estimated dark matter halo 
mass of the galaxies are Mdm = 1.2 x 10^^ Mq and 
Mbm = 1.0 X 10" Mq, respectively (Section [121) • Defin- 
ing the baryon mass fraction as /b = Mbtot/(-^^DM + 
Mb tot), we obtain /b w 0.11 for NGC1961 and /b « 0.09 
for NGC6753. 

The cosmic baryon mass fraction has been pre- 
cisely measured by the Wilkinson Microwave Anisotropy 
Probe (WMAP), yie lding A wmap = 0.171 ± 0.009 
(jPunklev et al.ll2009l ). This value is significantly larger 
than that observed in NGC1961 and NGC6753, imply- 
ing that '--^30 — 50% of the baryons are missing. This 
result is in agreement with obse rvations of othe r disk 
galaxies, and galaxy groups (see iBregmaiil 120071 for a 
review). However, as a caveat, we emphasize again the 
large uncertainties associated with the determination of 
the hot gas mass within the virial radius. This plays a 
particularly important role given that ^50% of the total 
baryon mass is in the form of hot gas in the (yet un- 
detected) outer parts of the halo. In principle, the real 
gas density profile beyond ~0.15r2oo may significantly 
deviate from the extrapolated ones, yielding a very dif- 
ferent total hot gas mass. For example, if the hot gas 
mass of NGC1961 and NGC6753 is factor of about 2-3 
times more than that estimated based on the density 
profiles, the resulting /b would be approximately consis- 



tent with cosmic baryon mass fraction. To undoubtedly 
measure the baryon budget of these (and other) spiral 
galaxies, further deep X-ray observations are required, 
which probe their hot gas content out to larger radius. 

6.3. Iron bias and abundances 

The major parameters of the X-ray coronae, listed in 
Tables [3] and HI are deduced from the best-fit spectra. 
However, two well-known complications are associated 
with the analysis of th ermal X-ray spectr a . One of them 
is th e iron bias (e.g. iBuote et al.l l2000t iBogdan et al.l 
|2012| ). the other is the degeneracy between the emission 
mea sure and the abun dance of the thermal component 
(e.g. iDavid et al.l[200l . 

In principle, the observed metal abundances may 
partly be influenced by the iron-bias, which appears if 
either a multi-temperature plasma or a temperature gra- 
dient is fit with a single temperature model. If the hot gas 
around NGC1961 and NGC6753 has a complex thermal 
structure, the obtained metal abundances may be under- 
estimated. The degeneracy between the emission mea- 
sure and the abundances of the thermal model may have 
an impact on the observed gas parameters. For example, 
if the abundance of the X-ray gas is not strictly sub-solar 
but is about Solar, i.e. ^10 times higher than measured, 
the gas densities, gas masses, and cooling times must be 
reduced by factor of about ~3.3. 

To address these points, the thermal and abun- 
dance structure of the hot gas must be better ex- 
plored. However, based on the presently available XMM- 
Newton data set, the application of more complex (e.g. 
multi-temperature) models is not feasible. Therefore, 
additional deep X-ray observations of NGC1961 and 
NGC6753 are needed to settle this issue. 

6.4. Future prospects 

Despite the extensive search in the last decades, we 
are the first to detect and characterize extended X-ray 
coronae beyond the optical disks of two normal spiral 
galaxies, NGC1961 and NGC6753. However, to further 
probe the X-ray coronae and constrain structure for- 
mation models, it is essential to identify X-ray coronae 
around other massive spiral galaxies. 

The large effective area of XMM-Newton makes it es- 
pecially suitable for the detection of hot gaseous coro- 
nae. Indeed, by combining the data of the three EPIC 
cameras, it is possible to collect a fairly large number 
of photons in moderately deep observations. Moreover, 
the large field-of-view allows to explore the coronae out 
to sufficiently large radii. A disadvantage of XMM- 
Newton is its relatively high instrumental background 
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level, which results in notable systematic uncertainties, 
thereby hindering the study of faint X-ray emission. For 
example, in NGC1961 and NGC6753 the systematic un- 
certainties dominate beyond ^^60 kpc radius, hence our 
understanding of their X-ray coronae would not signif- 
icantly improve at larger radii even with much deeper 
observations. Chandra has a lower level of instrumental 
background, therefore it is ideal to explore the coronae 
at large galactocentric distances. However, due to the 
smaller effective area of Chandra it takes significantly 
longer to collect a meaningful number of photons. 

The upcoming eROSITA r nission will perform a four 
year all sky-survey program (jPredehl et al.|[2"010[ ). dur- 
ing which a mean exposure time of 2.0 ks will be ob- 
served. This mission will play a major role in identifying 
the X-ray coronae of all nearby massive spiral galaxies. 
However, due to the relatively short exposure time of the 
eROSITA survey, follow-up observations will be required 
to characterize the X-ray coronae, which can be done 
with Chandra, XMM-Newton, or eROSITA itself. 

On a longer time scale, the Square Meter Arcsecond 
Resolution X-ray Telescope (SMART-X) mission will of- 
fer an outstanding possibility to characteri ze hot X-ray 
coron ae around massive spiral galaxies (jVikhlinin et al.l 
I2012f) . With its large field-of-view (22' x 22') and large ef- 
fective area (2.3 m^ at 1 keV) the SMART-X Active Pixel 
Sensor Imager will collect 5500 counts around NGC1961 
in a 10 ks observation in the 0.7 — 2 keV energy range 
from the region between 25 — 45 kpc. This will dras- 
tically increase the signal-to-noise ratio, thereby allow- 
ing to study the X-ray coronae in unprecedented details. 
Additionally, with the use of the X-ray Microcalorimeter 
Imaging Spectrometer aboard SMART-X the outer re- 
gions of the X-ray coronae can also be explored. Indeed, 
with the microcalorimeter the emission arising from the 
soft Galactic emission can be differentiated from the hot 
X-ray gas corona. Thus, it will be feasible to investi- 
gate the X-ray coronae of a large sample of disk galaxies, 
thereby further constraining galaxy formation scenarios. 

7. CONCLUSIONS 

In this paper, for the very first time, we detect and 
characterize the hot gaseous X-ray coronae around two 
"normal" spiral galaxies, NGC1961 and NGC6753, based 
on moderately deep XMM-Newton observations. Since 
luminous X-ray coronae are a fundamental prediction of 
structure formation models, we use our observations to 
probe hydrodynamical structure formation simulations. 
Our results can be summarized as follows. 

1. We detect gaseous X-ray coronae around NGC1961 
and NGC6753 out to ^60 kpc, which significantly 
exceeds the optical radii of the galaxies. The spa- 
tial distribution of the hot gas is fairly round and 
symmetric, indicating that the X-ray gas resides in 
hydrostatic equilibrium in the potential well of the 
galaxies. 



2. We find that the observed characteristics of 
the coronae are fairly similar in NGC1961 and 
NGC6753. In particular, within the (0.05 - 
0.15)r2oo region we measure gas temperatures of 
kT 0.6 keV and abundances of ^^0.12 Solar. The 
cooling time scale of the gas is ~50 Gyrs, imply- 
ing its quasi-static nature. The bolometric X-ray 
luminosity of the hot gas within the same region is 
~6 X 10""^ erg s"'^ for both galaxies. 

3. We derive the baryon mass fraction of NGC1961 
and NGC6753, and obtain /b w 0.11 and /b « 
0.09, respectively. These values are lower than 
the cosmic baryon mass fraction, indicating that 
~30 — 50% of the baryons are missing from these 
galaxies. However, large uncertainties are associ- 
ated with the derivation of the total hot gas mass 
within the virial radius, which could bias the ob- 
tained values of fh- 

4. The observed characteristics of the X-ray coronae 
are compared with the moving mesh code AREPO 
and the SPH-based code gadget. The imple- 
mented subresolution physics and the gravity solver 
are identical in the two codes, but they use dif- 
ferent methods to solve the equations of hydrody- 
namics. We find that, while neither model gives a 
perfect description, the AREPO code better repro- 
duces the luminosities, gas masses, and abundances 
observed around NGC1961 and NGC6753. Addi- 
tionally, the shape of the observed density profiles 
are also in better agreement with the moving mesh 
code within ~0.4r2oo- However, neither model in- 
corporates efficient feedback from AGN and super- 
nove, which could significantly alter the simulated 
properties of the coronae. 
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